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Models of eternal inflation predict a stochastic self-similar geometry of the universe 
at very large scales and allow existence of points that never thermalize. I explore the 
CN . fractal geometry of the resulting spacetime, using coordinate-independent quantities. 
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' The formalism of stochastic inflation can be used to obtain the fractal dimension of 

(N 

' the set of eternally inflating points (the "eternal fractal"). I also derive a nonlinear 
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branching diffusion equation describing global properties of the eternal set and the 
probability to realize eternal inflation. I show gauge invariance of the condition 

(N 

' for presence of eternal inflation. Finally, I consider the question of whether all 

oo ; 

■ thermalized regions merge into one connected domain. Fractal dimension of the 

eternal set provides a (weak) sufficient condition for merging. 
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I. INTRODUCTION 

In models of cosmological inflation, amplification of quantum fluctuations of the scalar 
field (the inflaton) may result in eternal self-reproduction of inflating domains, or eternal 
inflation 0, ^ In this case, stochastic evolution produces a spacetime in which ther- 
malization never happens globally. At arbitrarily late times there exist vast regions of space 
where inflation still continues; however, along any given comoving world-line, thermalization 
is certain to be reached (with probability 1) at a sufficiently late proper time. This some- 
what paradoxical situation is possible when self-reproduction of inflating regions proceeds 
faster than the drift of the inflaton field toward thermalization within these regions. 

To study the spacetime resulting from eternal inflation, one may use a stochastic formal- 
ism [0, 1^, 1^, 1^ that describes the probability distribution of the inflaton field 0, coarse-grained 
in horizon-size volumes of space. In simplest one- field models of inflation, the time- dependent 
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distribution of physical volume Py {(f), t) satisfies a Fokker-Planck (FP) equation 

dPv 92 , d 



at a^=(«^"'-a^(''^-' + 3^^- « 

This equation^ is supplemented by appropriate initial and boundary conditions; D (0) and 
V {(f)) are diffusion and drift coefficients and H {(f>) is the effective expansion rate. The kinetic 
coefficients are related to the infiaton potential V (0) by 



i7rV{(P) ^ I dH 



ff(^) = ^^. !, = _.„=-__. (2) 

At late times, solutions of Eq. (||) are dominated by the eigenfunction ifjy {(f)) of the diffusion 
operator with the largest eigenvalue 7^/, 

Pv(0,t)~e^-Vy(0)- (3) 

If 7v' > 0, the total volume of the inflating domain / Py (0, t) d(f> grows exponentially with 
time, which indicates presence of eternal inflation. 

One problem with the present form of the stochastic formalism is its dependence on 
the choice of equal-time surfaces. By construction, the volume distribution at constant 
time Py (0, t) is not a generally covariant quantity, and it is no surprise that solutions and 
eigenvalues of the FP equation depend non-trivially on the choice of the time variable t. 
The resulting spacetime is highly inhomogeneous on large scales, and different choices of 
time slicing introduce significant biases into the volume distribution. For this reason, it has 
been difficult to draw unambiguous conclusions about eternally inflating spacetimes using 
stochastic formalism P, 0, If. One could even question the meaning of the condition for 
presence of eternal inflation, •yy > 0, since the eigenvalue •jy is itself gauge-dependent. 

Eternal self-reproduction of inflating regions gives rise to fractal structure of eternally 
inflating spacetime on very large scales. Previously, the fractal structure of the inflating 
domain has been investigated using distributions of on equal-time surfaces J^, |lT 



However, the fractal dimension dfr of the inflating domain obtained in this way is deter- 
mined by an eigenvalue of the FP equation and depends, in general, on the choice of time 
parametrization (except for the case when H {(f)) = const). 

Another open issue concerns the global topology of the eternally self-reproducing universe: 
do thermalized domains remain forever separated from each other by an inflating sea, or do 



^ Here and below we use the Ito factor ordering of the diffusion operator 
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all thermalized domains eventually merge into one connected domain surrounding inflating 
islands? If merging occurs, all other thermalized domains in the universe will eventually enter 
our horizon; otherwise, we shall remain forever causally separated from other thermalized 
domains by eternally inflating walls. An early work by Guth and Weinberg [0 showed that 
in the old inflationary scenario there is no percolation of thermalized bubbles. The same 
occurs in certain models of topological inflation |T^ with inflating domain walls. Numerical 
simulations of Linde et al. 0, Q suggest that thermalized regions merge in some chaotic 
inflationary scenarios. However, a conclusive general description of the topology of the 
thermalized domain is lacking. 

In the present work, I consider models of inflation that can be described by Eq. (|l]) or 
its generalizations for multiple scalar flelds. I demonstrate that certain global properties 
of eternally inflating spacetimes can be determined in a coordinate-independent way. The 
principal object under consideration is the "eternal fractal set" the set of all comoving 
points that never thermalize. The set E is deflned (Sec. |I|) independently of the choice of 
coordinates. I show (Sec. |1I D|) that its fractal dimension, dimi?, is a gauge-invariant quantity 
which coincides with the fractal dimension df,. of the inflating domain computed using the 
scale factor time variable. In Sec. |T| I study conditions for presence of eternal inflation. 
First, I show that the basic criterion, 7y > 0, is in fact gauge-invariant. Then I deflne 
the probability X (0) to have eternal points in a given comoving domain; this probability 
plays the role of density of points of E. The quantity X (0) is gauge-independent and 



satisfles a nonlinear branching diffusion equation derived in Sec. |III B| . This new equation 



is shown (Sec. pi D]) to provide a (sufficient) condition for presence of eternal inflation that 



is consistent with the flrst criterion 7y > 0. Finally, in Sec. |V| I consider the question of 
whether all thermalized domains ultimately merge and become connected to each other at 
sufficiently late times. I show that the probability for merging is constrained by the nonlinear 



diffusion equation of Sec. [Ill B| . A (weak) sufficient condition for merging, dimE < 2, is 
provided by the fractal dimension of the set E. I conjecture that the fractal dimension may 
provide a stronger criterion of merging of thermalized domains. I also generalize the results 
to a suitable class of global topological properties of the set E. 



4 



II. THE ETERNAL FRACTAL 

In the current picture of eternal inflation, certain regions of spacetime may expand by 
an arbitrary factor and, in particular, arbitrarily small sub-Planckian scales may become 
macroscopically large. There is evidence that spacetime cannot be considered classical below 



the Planck scale (see Ref. [15] for a review). Recently, it was shown that a hypothetical mod- 
ification of physics on sub-Planckian scales may change the power spectrum of inflationary 
perturbations and create a significant backreaction affecting the evolution [|r7|. In that 
case, the kinetic coefficients of Eq. (|I]) would have to be modified. However, the qualitative 
features of an eternally self-reproducing spacetime would remain the same. In the present 
work we shall assume that no new physics emerges from previously sub-Planckian scales in 
the course of inflationary self-reproduction. 

Let us consider the evolution of an arbitrary finite (comoving) volume of space chosen at 
some initial time, in a cosmological model with eternal inflation. Here and below we assume 
that the initial comoving volume at t = had a certain value 0o of the inflaton field and 
was physically of one horizon size (0o)- 

There is a non-vanishing, if small, probability that the entire comoving region will ther- 
malize at a certain finite time, because fluctuations of the scalar field could have acciden- 
tally cooperated to drive (f) toward the end point of inflation 0^, everywhere in the region. 
One could say in that case that eternal inflation has not been realized in that region.^ If, on 
the other hand, eternal inflation has been realized, then for an arbitrarily late time t there 
will be regions that are still inflating at that time. 

In that case, there must exist some comoving world- lines that never enter a thermalized 
domain. To demonstrate this, one could choose a monotonically increasing sequence of time 
instances, t„ = nH~^, n = 0,1, 2, for each n there must exist a point x„ surrounded by 
a horizon-sized domain that will be still inflating at time t^. (As t„ grows, these domains 
become progressively smaller in comoving coordinates, due to expansion of space.) An 
infinite sequence of points Xn on a finite (compact) comoving volume of 3-dimensional space 
must have accumulation points, i.e. there must be at least one point x* such that any 
arbitrarily small (comoving) neighborhood of contains infinitely many points Xn from 



^ Here we do not consider the possibility of spontaneous return of thermalized regions to inflation in "recy- 
cling universe" scenarios [ISl. 
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the sequence. It is clear that the comoving world-hne at an accumulation point x^, cannot 
reach thermalization at any finite time: if it did, there would exist a comoving neighborhood 
around which thermalized at that time, and this contradicts the construction of the point 
x^,. We shall refer to points x* that never thermalize as "eternal points" and define the set 
E of all such points, as a set drawn on the spatial section of the comoving volume at initial 
time. (One could also imagine an infinitely dense grid of comoving world-lines starting at the 
initial surface, with the set E consisting of all world-lines that never reach thermalization.) 
The set E for a comoving region is not empty if eternal infiation has been realized in that 
region. 

Different histories and different choices of the initial comoving volume will generate differ- 
ent sets E. Since E is a stochastically generated set, one may characterize it probabilistically, 
by finding probabilities for the set E to have certain properties. The underlying stochastic 
process is the random walk of the infiaton and, since it is a stationary stochastic process 
that does not explicitly depend on time, probabilities of any properties of the set E depend 
only on (pQ. For instance, below we shall denote by X (0o) the probability for the set E to 
be non-empty if the initial comoving volume has horizon size and has a given value = 0o 
of the infiaton field. 

During stochastic evolution, there will be (infinitely many) times when the infiaton re- 
turns to the value 0o iii some horizon-sized domain; each time the distribution of properties 
of the subset of E within those domains will be the same as that for the whole set E. In 
this way, the set E naturally acquires fractal structure. 



The comoving volume P (0, t) of infiating regions with field value at time t is described 
by the FP equation 



A. 



The fractal dimension 




(4) 



with the kinetic coefficients from Eq. (^. Its late-time asymptotic solution is 



(5) 



The total comoving volume of infiating regions / P (0, t) dip decays exponentially rapidly at 
late times, and any pre-selected comoving worldline eventually reaches thermalization with 
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probability 1. In other words, the comoving volume occupied by eternal points vanishes 
and the eternal set E has measure zero, as a subset of 3-dimensional (comoving) space. Its 
fractal dimension must therefore be less than three. 

The fractal set E can be thought of as an infinite-time limit of the (fractal) inflating 
domain considered in Refs. |10|, Since the volume of E vanishes, we need a slightly 
different definition of fractal dimension than that of Ref. |TD|. The definition of fractal 
dimension we would like to use is the so-called "box fractal dimension" , which in our case is 
equivalent to the Hausdorff-Besicovitch fractal dimension . The "box fractal dimension" 
is defined for a non-empty subset of a finite region of a (i-dimensional space. Rectangular 
coordinates are chosen in the region to divide it into a regular lattice of d-dimensional cubes 
with side length /. (The cubical shape of the lattice is not essential to the definition of the 
"box fractal dimension".) Then, the "coarse-grained volume" V (l) of the set 5* is defined 
as the total volume of all cubes of the lattice that contain at least one point from the set S. 
If the set S has fractal dimension less than d, one expects V (/) to decay at / — > as some 
power of I. The fractal dimension is defined by 

dims = d-\im^^^^^. (6) 

For example, if the set 5* is made of a finite number of line segments, then the coarse-grained 
volume at scale / will be equal to the volume of narrow tubes of width / surrounding the 
line segments. Then V (/) ~ l'^'^ and the fractal dimension defined by Eq. (|^) would come 
out to be 1. 

Below we shall see that if the eternal set E is not empty, it has a well-defined fractal 
dimension dimii^ which is independent of the choice of spacetime coordinates, and we shall 
give a procedure to compute dimii^ in a given inflationary model. 



B. Topology of the eternal fractal 

Some simple topological properties of the eternal set E follow directly from its construc- 
tion. 

The set E is topologically closed: if a point x is not an eternal point, x^E, then there 
exists a neighborhood of x which also does not belong to E. This is so because the point x 
must have thermalized at a finite time t together with a horizon-size region around it. 
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The set E does not contain isolated points: in any neighborhood of a point x^E there 
exists, with probabihty 1, another point from E. To show this, suppose that a neighborhood 
of a point xeE contains no other points of E. Define to to be the time when this neighbor- 
hood grows to horizon (or larger) size and let Hq be the Hubble constant at point x at that 
time. Consider an infinite sequence of horizon-sized inflating regions i?„ surrounding the 
point X at times tn = + ''T'Hq^, n = 0,1, 2, by assumption, the evolution of the region 
Rn at each e-folding is such that Rn expands into inflating subdomains, of which one is to 
become -R^+i and all others contain no eternal points. There is a certain probability p < 1 
for a horizon-sized region i?„ containing one eternal point to have such an evolution during 
one e-folding; but the probability for all regions n = 0,1, 2... to have the same evolution 
is the product of infinitely many factors p, which is zero. 

By extension, any neighborhood of x G -E contains, with probability 1, infinitely many 
other points from E. 

Since E has measure zero, any given comoving point has a priori zero probability to 
belong to E. Therefore, any fixed line segment or surface has zero probability to consist 
entirely of points from E. 



C. Gauge invar iance of fractal dimension 



Previous calculations 0, [T^) [HI concerned the fractal dimension of the infiating domain 
at a finite time and assumed a certain time slicing (gauge). The results depended on the 
choice of the time coordinate (except for the case of Ref. |jll| where expansion of infiating 
domains was homogeneous and all choices of time coordinate are equivalent). Here, instead 
of the set of points that are infiating at a given time, we consider the set E of eternal points 
as a subset of comoving spatial section at initial time. Now we shall show that the fractal 
dimension of the set E is invariant under smooth changes of spacetime coordinates. 

The set E can be thought of as a locus of eternal points drawn in comoving spatial 
coordinates on the initial space- like slice (at t = 0). The construction of the set E requires 
to know whether a given point has thermalized at any finite time t > 0. Therefore the 
eternal set does not depend at all on the choice of time slicing of the spacetime to the future 
of the original slice. Its properties may, at most, be functions of the initial value of at 
t = 0. 
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Further, one can show that the fractal dimension of a set drawn in a finite region of 
space remains unchanged under any smooth change of coordinates in that region. First, 
the fractal dimension of a set is unchanged under a linear transformation of coordinates 
with a non-degenerate Jacobian J, because the coarse-grained volume V (/) will simply gain 
a factor J and this will not change the fractal dimension defined by Eq. (^. Given an 
arbitrary smooth coordinate transformation, we can choose a sufficiently small scale Iq at 
which the coordinates are locally changed by a constant linear transformation. We can 
divide the initial region into subdomains of size li < Iq and find V (/) for Z < /i by adding 
the volumes Vi (/) calculated over these subdomains. Under the coordinate change, each 
volume Vi (/) will be multiplied by a corresponding Jacobian Jj, but the Jacobian in a given 
subdomain does not depend on / for sufficiently small h. Therefore, in the limit / the 
dominant dependence of V (/) on I of the form V (Z) ~ P will be unchanged under coordinate 
transformations. 

Finally, the fractal dimension of the eternal set E is the same when computed in an 
arbitrarily small neighborhood of any eternal point. This is because the set E has no 
isolated points and any point x & E is surrounded by infinitely many other eternal points 
lying arbitrarily close to x. The statement can be justified more formally as follows. In the 
next subsection we show (without using the present statement) that the fractal dimension 
of E is, in fact, independent of the initial value of in the initial horizon-sized region. If 
we choose a small comoving neighborhood U {x) of an eternal point x, there will be a time 
ti when U (x) grows to horizon size. Since U (x) contains an eternal point x, it could not 
have thermalized by the time ti. Therefore at time ti the neighborhood U {x) is undergoing 
infiation and it can be taken as the initial comoving region with a certain value of 0. The 
future evolution of any infiating horizon-sized region is statistically the same and therefore 
the fractal dimension of the subset of E computed in the neighborhood U (x) starting at 
t = ti is the same as that of the full set E. 

Unlike other fractals occurring in nature, the fractal set E has no short-distance cutoff: 
it retains its fractal structure at all scales (in comoving space). 



9 



D. Fractal dimension from Fokker-Planck equation 

In this section we show that the fractal dimension of the eternal set E can be obtained 
from the comoving-volume probability distribution of the field in a special gauge. This 
probability distribution is a solution of the comoving-volume FP equation in that gauge. 
The dominant eigenvalue of this FP equation determines the fractal dimension of the set 
E and coincides with the dominant eigenvalue of the FP equation in the scale factor time 
variable. 

We introduce a special time variable 6 related to the scale factor a by 

g(t,x) = lna(t,x)+ln ^^/;''\ (7) 

-no 

The meaning of this variable is the physical expansion scale a {t, x) relative to the local 
horizon scale H~^; a surface of equal 6 consists of spacetime points (t,x) where the initial 
horizon size Hq^ had expanded to times the current horizon size at the same point. One 
may call 6 a "horizon scale factor" time variable. 

The time variable 6 is not strictly monotonic: 9 = a/a (dots denote proper time deriva- 
tives) and if H is decreased rapidly enough so that the accelerated expansion locally stops, 
the "time" 6 will decrease despite growth of the scale factor a. However, this is expected 
to happen only near the end of inflation, whereas in other regimes such large fluctuations 
of H are extremely rare and expansion is almost always accelerated. We shall assume that 
the time variable 9 is well-behaved throughout the allowed interval of and truncate that 
interval near the end of inflation if necessary. 

Below we shall see that in the gauge 9, the comoving volume probability distribution 
P*^^^ (6', 0) at late times 9 has the form 

pW(^^,0)^^W(0)e-^^ (8) 

where (—7) is the largest exponent dominating the solution at late times and ■0*^^^ is an 
appropriate eigenfunction. Assuming that Eq. (|]) gives the asymptotic solution at late 
times, one finds the fractal dimension of the set E to be 

dimE = 3-7. (9) 

This result can be derived from Eq. (P) as follows. Let = 3 be the dimension of 
space. An initial comoving region of horizon size with a constant field value 0o has linear 
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extent Hq^, where Hq = H ((^q), and we can choose comoving coordinates in the region to 
divide it into identical cubical subdomains of equal initial linear size IHq^, with / <^ 1. 
At a sufficiently late time 6i (/) when these comoving subdomains have grown to horizon 
size, there will be on average l^'^P^^^ {6i,(j))d(j) subdomains still inflating with field value 
within the interval [0, (f) + dcj)] . (We assume that eternal inflation has been realized in the 
initial region, otherwise the number of inflating subdomains would be zero at sufficiently late 
times.) The probability for a horizon-size region with field value to contain at least one 
eternal point is a function of (f) which we denote X (0). It follows that at time di there will 
be, on average, Z^'^/ P*-^-* {9i, (p) X {(p) d(f) subdomains that contain an eternal point. Since 
g0i _ ^-1 g^j^^ 0^ jg lo^j-gQ ^Qj- sjxiall I, we can use Eq. (§) to obtain the total (comoving) volume 
of those subdomains, 

V (Z) = e-^^i J (0) X (0) d(j) = const ■ P. (10) 

This is the coarse-grained volume of the eternal set at scale /. By definition of fractal 
dimension, dimE = (i — 7. 

The constant in Eq. (^) is nonzero if X (0) ^ and if the set E is not empty; otherwise 
we cannot assume that there are many infiating subdomains and estimate V (l) by its average 
value, as we have done in Eq. (|10|) . This constant also absorbs all dependence on the initial 
value of in the initial comoving region. It follows that the fractal dimension is independent 
of the initial value of and is well-defined as long as eternal infiation is realized in the domain 
of consideration. 

Now we need to determine 7. Instead of using a Fokker-Planck equation in the gauge 
6, we will show that 7 is equal to the dominant eigenvalue of the Fokker-Planck equation 
in the scale factor time variable. This will give a computational prescription to obtain the 
fractal dimension of the set E in a. given model of infiation. 

The distribution of volume P^^^ {9, 0) in the time gauge 9 is related to the distribution 
P (t'^^'^ , 0^ in the scale factor time gauge t^^-* = In a by 

P 0) d<f) = P(^) (^t(o) + In 0^ d<p. (11) 

To justify this, consider the combined comoving volume of all domains in which the value 
of the field is between 0i and 0i + dcp and scale factor is exp[t*^*^)]. This volume is 
P{t^^\(j)i)d(j). [The distribution P(t^^\(f)) as obtained from the FP equation is not nor- 
malized since / P{t^^\(j))d(j) < 1 is the total comoving volume of infiating regions in units 



11 



of the initial comoving volume.] Throughout all volume characterized by 0i and t^^\ the 
"horizon scale factor" 6 has the value e{t^°\ 0i) = t(°) + ln[if(0i)/i7o] + O {dcj)). Also, any 
volume in which = 9{t^^\ 0i) and G 0i + d(f)] must have scale factor exp[i(:*^°''] +0 (dcf)). 
Therefore, P{t^^\(pi)d(p is, up to O {d(f)^), equal to the total volume of all domains in which 
e = 0i) and G [01, 01 + d(f)]. This justifies Eq. (|Tl]). 

The same justification can be given for an analogous relation between physical-volume 
distributions in gauges 9 and t*-"-* , since the argument does not rely on the comoving nature 
of the volume to which the distributions relate. 

The FP equation for the comoving volume distribution in the scale factor time gauge t'-"-' 

is 

dm ~ 902 ^) 90 ^ ) ■ ^ ^ 

Here 

/)(o) = :^ ^(0)^ (13) 

87r2' 47rif rf0 ^ ' 

are the kinetic coefficients in this gauge. At late times t^^\ the solution of Eq. ([T^ ) is 

P(tW,0) ^e-^*'"V(0), (14) 

where 7 is the dominant eigenvalue of Eq. ([T2|). Late times ^ correspond also to large 
values of since H/Hq is bounded. It follows from Eq. ( pT]) that the distribution P^^^ has 
the asymptotic form 

p(^)(^^,0)^e--^(^) '^(0), (15) 

Therefore, the dominant eigenvalue 7 of the FP equation in the scale factor time gauge 
is the relevant quantity for the fractal dimension in Eq. (^. 

This result shows that the fractal dimension of the set E given by Eq. is the same as 
the fractal dimension df^ of inflating domain [l^, |lT| if the latter is computed on surfaces 



of equal scale factor. Although 7 is an eigenvalue of the FP equation in a particular gauge, 
the quantity d — 7 has a physical interpretation as the fractal dimension of a set defined 
independently of spacetime coordinates; eigenvalues of the FP equation in other gauges do 
not seem to have a gauge-invariant interpretation. 
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III. CONDITIONS FOR ETERNAL INFLATION 



We have distinguished two issues regarding the possibihty of eternal inflation: first, eter- 
nal inflation may be entirely disallowed in certain models; and second, in models where 
eternal inflation is possible, it may accidentally not have realized in certain domains. In this 
section, we first examine conditions on an inflationary model which make eternal inflation 
possible, and then we shall obtain the probability for eternal inflation to be realized in a 
given comoving region. The condition that the latter probability does not vanish is another 
criterion for presence of eternal inflation. Both criteria are shown to be independent of the 
choice of spacetime coordinates. Finally, we demonstrate that the two criteria for presence 
of eternal inflation are consistent with each other. 



A. Presence of eternal inflation 



The basic condition for possibility of eternal inflation in an inflationary model is 7y > 0, 
where 7y is the dominant eigenvalue of Eq. (0). However, this equation itself and the value 
of 7y are gauge-dependent. We shall now show that 7y > is in fact a gauge-invariant 
condition. 

The key idea of the proof is that if 7y = in some gauge it must also be equal to zero in 
all other gauges, and therefore 7y cannot be positive in one gauge and negative in another. 

Suppose that there exist two time variables, to and ti such that the dominant eigenvalue 
jv has opposite sign in the gauges defined by and ti. Introduce a one-parametric family 
of time variables to,, < a < 1 to interpolate between these time variables. 



dta _ ( dti 



(16) 



The new time variables ta are defined by integrating Eq. ([161) along comoving world-lines. 
The dominant eigenvalue 7y of the FP equation in gauge is a function of a. This function 
■jv (ct) is continuous, because the differential operator in the FP equation is self-adjoint, its 
dominant eigenvalue jv is non-degenerate and under a small change of a this eigenvalue 
must change by a small amount, found from standard perturbation theory. By continuity 
it follows that there exists a value ao such that 7y (ao) = 0. The FP equation in this time 
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parametrization has a stationary solution Py\ 

;|K^'')-|M"V3^^«4°'^0. (17) 

Here Da, Va and Ha are kinetic coefficients in the gauge ta- They differ from Dq {(f)), vq (0) 
and Hq {(p) in the gauge to by the factor (dta/dto)'^ . The factor dta/dto can be absorbed 
into Py\ and we obtain a stationary solution of the FP equation in all other gauges ta- 
Existence of this solution with eigenvalue in all gauges ta contradicts our assumption that 
the largest eigenvalue of the FP equation is negative in one of the gauges to or ti. 

We have found that the dominant eigenvalue jv of Eq. must have the same sign in 
all gauges and, therefore, presence or absence of eternal inflation can be judged by the sign 
of 'jv in any gauge. 

B. A diffusion equation for X {(p) 

Assuming now that 7y > 0, we consider the probability for eternal inflation to be realized 



in a given region. In Sec. [II D| , we have denoted by X (0) the probability for an initial inflating 
region of one horizon size, in which the scalar field has value (p, to contain at least one eternal 
point. It will be convenient to work with the complementary probability X (0) = 1 — X (0). 
In this section we find that X (0) satisfies a nonlinear branching diffusion equation [Eq. (p3D] 
that resembles the backward FP equation for the physical volume probability distribution. 
We derive this equation and show its independence of time parametrization. 

Here we shall use proper time t as the time variable. The random walk of the infiaton 
field is described by a Langevin equation which can be written as a difference equation 



(t + At) = (t) + V (0) At + ^^2D (0) At, (18) 

where ^ = ^ (t, x) is a normalized random "noise" representing fluctuations, and v (0), D (0) 
are the kinetic coefficients given by Eq. @. The evolution of in a horizon-size domain with 
field value = 0o at time t = gives rise to independent random walks of in each of the 
"daughter" horizon-size subdomains that were formed out of the original domain. After one 
e-folding, i.e. after time At = H~^, there are N = exp [3HAt) ^ 20 daughter subdomains of 
approximately horizon size. The stochastic process corresponding to ^ assigns a probability 
density P {^i, ■■.,^n) d^i...d^N for various sets of values of ^ in the daughter subdomains. 



14 



The probability X (0o) to have no eternal points in the region is equal to the probability 
to have no eternal points in any of its daughter subdomains. The evolution in inflating 
daughter subdomains proceeds independently, so the probability to have no eternal points 
in each of them is described by the same function X (0) evaluated at (t + At). (Although 
the daughter subdomains have slightly varying H (0) and are not exactly of horizon size, 
the correction due to this is negligible.) This gives an integral equation. 



X{4>o) = J d^^...d^NP{^l,...,^N) 



X n ^ (00 + (0o) At + ^i^2D (0o)At) . (19) 
1=1 ^ ^ 

Here we use the Ito interpretation of the Langevin equation, where v and D are evaluated 

at the initial point of the step, = 0o- 



Since correlations between different daughter subdomains are small |21|], we can approx- 
imate the probability density of by a product of independent identical distributions, 
P (^1, ■■.,^n) = P (6) ■■■P (^n)- Then Eq. (0) becomes 



X(0o) = 

/ diP (0 X (00 + vAt + iV2DKt) . (20) 



Rather than trying to solve Eq. (|20|) , we consider its limit for small values of HAt. Although 
the Langevin description of the random walk given by Eq. (|l^) is valid only for times 
if At > 1, we shall formally consider it to be valid at all values of At and take a partial 
derivative d/dAt of Eq. (^) at At = 0. The same formal limit At — is used to derive 
the FP equations such as Eqs. (|l]), (|12D in the standard stochastic formalism of inflation, 
and the same limits of validity apply to the new diffusion equation [Eq. (|23|) ] that will result 
from the present argument. 
We find 

^^ = = 3HXlnX 
dAt 

r dX (00 + t;At + eV2DAt) 

- 

xf-.|^). (21) 



2VAt / 
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We now expand dX / dip in Taylor series around = 0o and use the fact that ^ is a normahzed 
random variable, 

jd^p{o = h (0 = 0, {e) = h (22) 



to obtain an equation for X, 



d'^X dX 

D^ + v^ + 3HX In X = 0. (23) 



This equation needs to be supplemented with suitable boundary conditions at end of 
inflation and/or at Planck boundaries. The probability X to have no eternal points for a 
region which is near the end of inflation = 0^, should be 1. Domains reaching Planck 
boundaries (ppi will never thermalize and we could set X = at those boundaries; this is 
consistent with the viewpoint that super-Planck domains effectively disappear |^. Physically 
meaningful solutions of Eq. (|23| ) should vary between and 1. Therefore, the additional 
requirements are 

0<X<1; X(0,) = 1; X(0„) = O. (24) 

It is straightforward to verify that Eq. (pSj), unlike Eq. (|1|), is a gauge-invariant equation. 
A change of time variable, t t', with dt' /dt = T {t, 0) , will divide the functions D (0) and 
V {(f)), as well as the factor 3H (0) in the "growth" term, by T {t, 0). This extra factor T (t, 0) 
cancels and Eq. (^) remains unchanged. This is to be expected since the probability X (0) 
is defined in a gauge-invariant way and should be described by a gauge-invariant equation. 

Equation ( p3|) is a branching diffusion equation similar to the backward FP equation for 
the physical volume probability distribution in the Ito factor ordering, except for the extra 
factor InX that makes it nonlinear. An analogous diffusion equation can be derived for 
models of inflation with multiple scalar fields 0^, 

f) r)X F)X 

D^^ + Vk^ + ?,HX\nX = Q, (25) 



where 

1 dH 



are the new drift coefficients (in the proper time gauge). Boundary conditions will be 
X (0) = 1 along thermalization boundaries and X (0) = along Planck boundaries (if any). 

The qualitative behavior of the probability X (0) for a chaotic type inflationary model 
is illustrated in Fig. An approximate WKB solution in the fluctuation-dominated regime 
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FIG. 1: Probability X {(j)) to have no eternal points. 



can be obtained if we disregard the drift term v {(p) X' and write an ansatz 



X{<f)) = Ae~^^'^\ 



(27) 



where A is approximately constant and W{(j)) is a slow-changing function, \W"\ <^ {Wy 
that satisfies 



where (j)q is the boundary of fluctuation-dominated range of 0. It shows that the probability 
to have eternal points 1 — X (0) exponentially rapidly approaches 1 in the fluctuation- 
dominated regime. 

Since X > 0, there is always a certain (perhaps exceedingly small) probability for eternal 
inflation to not be realized within a given initial horizon-sized region. On the other hand, an 
inflationary model could be fine-tuned to entirely disallow eternal inflation, so that Eqs. (^3])- 
(0) have no solutions other than X = 1. Note that, in the absence of Planck boundaries, 
X = 1 is always a solution of Eqs. (p^)-(p^). This solution would mean that there are no 
eternal points anywhere. Equations (p3D-(p4D alone do not provide enough information to 
choose between the constant solution X = 1 and a nontrivial solution X (0). However, we 
know from considerations in Sec. |I| that, if eternal inflation is allowed in a model, then 
eternal points are possible, i. e. the probability to have eternal points is nonzero. Therefore, 
we should choose the solution X = 1 only when no other solution of Eq. (^31) that satisfies 
Eq. ( ^41) can be found. 



D{W'f -3HW = 0. 



(28) 



The solution is 




(29) 
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FIG. 2: Illustrative potential for eternal inflation. 
C. An exact solution 

For illustrative purposes, we solve Eq. (p3|) explicitly for a flat potential connected to 
two slow-roll slopes with negligible diffusion (Fig. |^). This potential has been considered in 
Ref. |]l9l as a qualitative illustration of eternal inflation. The allowed range of (p is divided 



into two regimes: in the fluctuation-dominated regime, 



< < 



J, the potential is flat 



and we take H (0) = Hq, D (0) = Dq and v (0) = 0. In the second regime (0]; < < 



or 0^ < 



< 



we take D (0) = and v (0) 7^ 0, to represent pure deterministic motion 



without fluctuations. 

We expect that X = 1 in the second regime, since all regions thermalize within finite 
time. In this regime Eq. ( ]23|) becomes 



vX' + 2HX\nX = 0. 
The general solution of this equation is 

X (0) = exp 

where 



(30) 



C 



a (0) = exp 



L «(0)' 

Hd(j)' 

, IF 



-Att 



(31) 



(32) 



is the scale factor along a slow roll trajectory, and C is an integration constant. The boundary 
condition X {(p^) = 1 forces C = and therefore X = 1 throughout the deterministic region, 
as expected. 

Turning now to the fluctuation-dominated regime, we need to solve 



DoX" + 3HoX\nX = 



(33) 



18 



with boundary conditions 

mi') = 1- (34) 

For simplicity, we take 0^ = 0^ = — 0^. Equation ([33| ) describes a Hamiltonian system of a 
particle with mass Dq moving in "time" in a potential 

UiX) = ^X^{2\nX-l). (35) 

This potential monotonically decreases from f/ = at X = to its minimum U = —3Ho/4 
at X = 1 (see Fig. Boundary conditions correspond to a "trajectory" X (0) starting and 
ending at X = 1 such that the total "travel time" is 20g. Such a trajectory is unique and is 
characterized by the lowest reached value Xq (the turning point) such that 

f"! dX\/Dn 



^0 ^2t/(Xo) - 2U{X) 



(36) 



The value Xo is the lowest probability to have no eternal points and is naturally achieved 
at = 0, in the middle of the fluctuation-dominated range of 0. In this model, Xq is a 
function of {(pq/Ho) which can be obtained numerically. For 0q ^ Hq, one can approximate 
U{Xq) ^ and obtain 

1 



Xq ~ 2 exp 



— - s 
2 



s ^ 2nV3^ > 1. (37) 
-no 



Numerical evaluation shows that this asymptotic expression holds to 5% accuracy for 
2TT(j)g/Ho > 2 which corresponds to Xq < 10^'*. As expected, the probability X = 1 — X 
to have eternal points is very close to 1 for regions that start in the fluctuation-dominated 
regime. 

The opposite case, Xq ~ 1, corresponds to motion near the minimum of the potential 



f/(X). Since the period of such motion is approximately independent of amplitude, Eq. (|36D 
cannot be satisfied unless 

0. > 0c = (38) 

For 0g < 0c, Eq. (|33|) has no solutions except X = 1. This indicates absence of eternal 
inflation when the fluctuation-dominated range of is sufficiently narrow. 

It is straightforward to check that for the inflaton potential V (0) of Fig. 0, Eq. (|l|) with 
boundary conditions Py {4>±qyt) = admits a non-negative eigenvalue ■jv if and only if 
Eq. ( p8D holds. In the next subsection we shall see that this result is not a coincidence. 
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FIG. 3: Effective potential of Eq. (BB) 



D. Consistency of conditions for eternal inflation 



We have given two gauge-invariant criteria for presence of eternal inflation: positivity of 
the dominant eigenvalue jv of Eq. (|l|) and existence of a non-trivial solution of Eqs. (p3D -(p4|). 
Since the relevant diffusion equations are somewhat dissimilar, a natural question is whether 
these two criteria are equivalent. A positive answer is suggested by the analytic example of 



the previous section. Here we shall see that Eqs. (23)-(24) cannot have a nontrivial solution 
X {(p) ^ 1 if the FP equation has a non-positive dominant eigenvalue.^ 

The idea of the proof is to consider a suitable non-negative function of 0, integrate it 
over (f) and obtain 7v > ClS db result. 

Let C be the differential operator of Eq. (|l|), 



CPv 



o2 f\ 

—r (DPv) - ^ (vPv) + 3HPv. 



(39) 



The spectrum of eigenvalues of C is bounded from above, and the eigenfunction ipy (0) of 
Eq. corresponding to the largest eigenvalue is everywhere positive, as a ground state 
of a self-adjoint operator Suppose that a non- trivial solution X (0) ^ 1 of Eqs. (|2|)-(|2^ 
exists; then the following integral must be positive. 



r{i-x 



ipv (0) d(j) > 0. 



(40) 



^ I do not yet have a proof of existence of a nontrivial solution X {(f)) when 7y > 0, which is needed to 
rigorously demonstrate equivalence of these two criteria. 
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Here 0i_2 are left and right thermalization boundaries, 02 > 0i- Since 

1 r<l>2 



h hi ^ ' 

we would obtain the desired inequality 7y > if we prove that 

h = J^' (l - X) Cijydcf) > 0. 
Integrating Eq. p^) by parts and using boundary conditions X (0i,2) = 1, we find 



(41) 



(42) 



1-X 



(43) 



Here C* is the operator conjugate to C. Using Eq. ([23|), we obtain 



1-X 

d^x dx 



-D 



V— + 3H 1-X 



3H (l- X + X\nX) > 0. 



(44) 



The integral / '0y£*[l — X]d(f) is positive because its integrand is nonnegative and not ev- 
erywhere zero. The boundary term DipyX' is non-negative at 02 and non-positive at 0i. 
This proves Eq. (^2]) and therefore '-jy must be positive. 

It follows that if the dominant eigenvalue 7^ is zero or negative, a nontrivial solution 
X (0) ^ 1 cannot exist. 



IV. MERGING OF THERMALIZED REGIONS 



In the previous section, we found a function X (0) that describes the probability to realize 
eternal inflation. Here we show that the same function gives bounds on probabilities of other 
global properties of the set E. In particular, we consider the issue of whether all thermalized 
domains ultimately merge and become connected to each other. 



A. Global properties of the eternal set 



Observe that the nonlinear diffusion equation for X (0) was derived using the argument 
that an inflating domain contains no eternal points if and only if none of its daughter 
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domains, after one e-folding, contain eternal points. It is clear that the same derivation 
would apply for the probability Y (0) to realize some other property Y of the set E, as 
long as that property holds for a given horizon-size inflating domain if and only if it holds 
for all daughter domains, and given that the property Y holds for a thermalized domain 
{Y (0^,) = 1). These conditions define a certain class G of "global" properties that are binary 
alternatives, functions of the subset of E inside a given domain. 

We find that the probability to realize any property from this class G is described by 
Eqs. (^3D-(p4|). The property for the set E to be empty, which occurs with probability 
belongs to this class; other examples would be the property for the set E to not 
contain any continuous line segments, or any continuous 2- dimensional surface segments. 
Let us denote by Sn the property of the eternal set to contain no fragments of continuous 
72- dimensional submanifolds, n = 0,1, 2, 3. 

The universal applicability of Eq. (p3D to all global properties of class G can be interpreted 
as a statement that probabilities for global properties of the fractal set E must be "fixed 
points" under the "renormalization" (rescaling) of the set, and all such fixed points must 
satisfy Eq. (|2|). 

In presence of eternal inflation, as we have seen, Eqs. p3|)-(|2^) admit two solutions, 
namely X = 1 and X = X {(p) ^ 1. Equations (^)-(^) alone do not give enough in- 
formation for us to select one of these two solutions. To find the probability distribution 

Y {(f)) for a particular property Y, other considerations specific to that property and to the 
given inflationary model are needed. However, we can conclude a priori that one of these 
two alternatives must be realized for any given property Y from the class G. This implies 
that a property from the class G either always holds, Y (0) = 1, or holds with probability 

Y {(f)) = X{(f>). (Here X {(f>) is a fixed function, independent of the property Y.) For in- 
stance, it would be enough to show that Y {(f)) > X (0) for some property Y to conclude 
that r (0) = 1. 

Since the same function X (0) describes the probability to have no eternal points, a 
natural question to ask is whether the occurrence of a given property in a comoving domain 
is correlated with presence of eternal points in the same domain. Consider a property Y of 
class G and suppose that we know that Y does not always hold; it follows that the property 

Y holds with the probability Y (0) = X (0) ^ 1 in a horizon-size domain with value of 
the inflaton field. If a horizon-size domain has no eternal points (property X), then the 
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property Y must hold in this domain. Therefore the probabihty for Y to hold if there are no 
eternal points is Prob{X CiY) = Prob{X) = Prob{Y), while the probability for Y to hold 
in presence of eternal points is Prob{X (lY) = Prob{Y) — Prob{X (lY) = 0. 

We find that, for any property Y of class G, one of the following two alternatives must 
hold (in a given inflationary model): either Y is always true for any horizon-size domain, or 
Y is true if and only if the domain contains no eternal points. In other words, a property Y 
of class G either always holds or is perfectly anti-correlated with presence of eternal points 
and, by extension, perfectly correlated with other properties of class G. For example, the 
property Sq holds in a domain only when the domain contains no eternal points, while S3 
always holds. In principle, details of the model should determine which of these alternatives 
is realized for a particular property Y. 



B. Bounds on merging of thermalized regions 

In this section, we give a bound on merging of thermalized regions based on the fractal 
dimension of the eternal set E. 

The set of all comoving points that will eventually thermalize is a complement of the set 
E. Thermalized regions will not merge if the set E encloses pockets of thermalized comoving 
space. To enclose a region of a 3-dimensional space, the set E must contain a continuous 
boundary of the enclosure which is at least a 2-dimensional surface, so the fractal dimension 
of E must be at least 2. If the fractal dimension dimE of the set E turns out to be less 
than 2, all thermalized domains must be merged. If, on the other hand, dimE > 2, then 
the question remains unresolved, because there exist sets of fractal dimension up to 3 which 
nevertheless do not enclose any interior domains. 

The condition dimii^ < 2 gives a topological bound on the possibility of merging of 
thermalized regions. However, this bound may be rather weak because in typical inflationary 
scenarios dimE = 8 — 7 with 7 <^ 1. For instance, in a potential with a fiat maximum at 
= 0, an estimate of Ref. gives 

^^rm«i. (45) 

^ 167r V (0) ^ ^ 

In this case the condition 7 > 1 is not satisfied and the question of merging remains open. 
One may distinguish merging of all thermalized regions into one connected domain from 
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percolation of thermalized regions, that is, formation of at least one infinitely large thermal- 
ized cluster (of infinite comoving volume). Percolation is a weaker condition, since merging 



entails percolation but not vice versa. Guth and Weinberg have considered the "old" 
scenario of inflation with bubble nucleation and gave bounds on the nucleation rate e per 
Hubble 4-volume for the bubbles (representing thermalized regions) to percolate. They have 
found that e < 1.1 x 10^^ is a sufficient condition for non-percolation of bubbles, while 
e > 0.24 guarantees percolation. In the model of bubble nucleation, the fractal dimension 



of the inflating domain at constant proper time is [H 



dfr = S- Ane/3. (46) 

Since in this model the expansion of inflating domains is homogeneous with H = const, time 
parametrization by scale factor is equivalent to proper time parametrization, and the fractal 
dimension of the set E has the same value, dimE = dfr- It is interesting to note that the 
sufficient condition given in Ref. for percolation of bubbles, e > 0.24, approximately 
coincides (to 2 decimal digits) with the topological bound dfj. < 2 which guarantees merg- 
ing.^ I conjecture that, in models of scalar field inflation, an upper bound on the merging 
probability may be established using the fractal dimension of the set E: namely, merging 
does not occur if dimE is sufficiently close to 3. 

The property that all regions that eventually thermalized (all points outside the set 
E) are merged into one connected region is not a property of class G, because merging 
of thermalized regions in a domain does not entail merging in all of its subdomains (it is 
possible that some subdomains are divided by fragments of eternally inflating walls which, 
nevertheless, do not globally enclose an interior region). We can consider a weaker property 
that belongs to the class G: the property 5*2 for a region to contain no eternally inflating 
2-dimensional surfaces. This property is sufficient (but not necessary) for merging. On the 
other hand, if a domain has a nonzero probability to contain fragments of 2-dimensional 
eternally inflating surfaces, it will contain, with the same probability, an inflnite number of 
arbitrarily small such fragments. If this probability has been realized, it is likely (but not 
certain) that merging of thermalized regions will be impossible in that domain. 

From considerations in the previous section we may conclude that, in a given inflationary 
model, either the property 5^2 always holds (and as a consequence all thermalized domains 



^ I am grateful to A. Vilenkin for bringing this point to my attention. 
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will always merge), or the property 5*2 never holds in regions where eternal points are present. 
This would provide another bound for the probability of merging if the probability to realize 
the property 5^2 were determined. 

I have been able to show that in any model of scalar field inflation described by Eq. ([1|), the 
eternal set E cannot contain any differentiable line segments or fragments of different iable 
surfaces (details of the proof will be given elsewhere). The property of E to not contain any 
differentiable line segments also belongs to the class G and its probability can be shown to 
exceed X (0); then it follows from considerations of the previous section that this property 
must hold for all regions. However, merging of thermalized domains does not follow from 
this result, because the fractal set E may well contain continuous but not differentiable lines 
or surfaces. 



V. SUMMARY AND DISCUSSION 

In this work I have explored the fractal geometry of an eternally inflating universe at very 
large scales, in the formalism of stochastic inflation. I defined a fractal set E consisting of 
all eternally inflating points in comoving space. The definition of set E is independent of the 
choice of spacetime coordinates. Using the Fokker-Planck equation for comoving volume, 
I found the fractal dimension of the set E. It coincides with the fractal dimension dfr of 



the inflating domain as defined in Ref . |jTO[ , if the latter were computed on surfaces of equal 
scale factor instead of equal proper time. In this way, the gauge-invariant construction of 
the fractal set E resolves the issue of gauge dependence that plagued previous calculations 
of the fractal dimension. 



In a recent work, Bousso p2| considered infinite self-reproduction (proliferation) of de 
Sitter space due to nucleation of black holes. This effect may also create a fractal structure 
of the spacetime on very large scales. Although the proliferation effect involves topology 
change which is absent in the usual picture of eternal inflation, the fractal structure of the 
resulting spacetime can be investigated using the methodology developed above. 

I have also examined the conditions for presence of eternal inflation in a given inflationary 
model. The original criterion was an unbounded growth of physical volume of inflating 
domains, which is equivalent to the condition 7y > imposed on the largest eigenvalue 7y 
of the FP equation for the physical volume [Eq. (|l])]. I have demonstrated that this condition 
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is gauge-invariant by proving that the sign of the dominant eigenvalue of Eq. ([^) must 
be the same in all gauges. Then I found the probability X (0) for an initial horizon-size 
region with scalar field value to contain at least one eternal point. This probability is a 
solution of a (gauge-invariant) nonlinear diffusion equation [Eq. (0)] derived in Sec. [Ill Bl 
and analyzed in the following sections. To show that the new nonlinear diffusion equation is 
consistent with the existing formalism, I checked that the probability X (0) vanishes when 
eternal inflation is not present, in agreement with the criterion 7y > 0. 

Finally, I have investigated the issue of merging of thermalized regions into one connected 
domain, as observed at arbitrarily late times. I obtained a topological bound on merging: 
all thermalized regions are guaranteed to merge if the fractal dimension of the set E is less 



than 2. This bound agrees with the previously found sufficient condition e > 0.24 |T2[ for 



percolation in the model of bubble nucleation (in that model, Eq. (46) relates the nucleation 
rate e to the fractal dimension of the set E). In Ref. [0, a sufficient condition for non- 
percolation, e < 1.1 X 10~^, was also found. This suggests a possibility of a stronger bound 
on merging, namely that merging does not occur if the fractal dimension is sufficiently close 
to 3. 

The stochastic formalism also constrains the probability of merging in a different way. 
Since the fractal structure of the eternally inflating set is preserved on all scales, the probabil- 
ity to realize any global property of the set must be invariant under rescaling. Equation ( pSf ) 
expresses the scale invariance of the probability X (0) for a horizon-size region to thermalize, 
and the same equation applies to probabilities to realize other global properties. I have de- 
fined a suitable class G of global properties described by Eq. (^); for example, the property 
5^2 for a domain to contain no eternally inflating fragments of 2-dimensional surfaces is a 
property of class G. The property 5*2 is a sufficient condition for merging of thermalized 
domains (which is itself not a property of class G). I found that the probability to realize 
a property of class G is either identically equal to 1 or is described by the function X (0). 
This provides a constraint on the probability of merging. However, Eq. (p3|) alone does not 
provide enough information to resolve this alternative. This is to be expected, because the 
FP equation concerns only the total volume of space with given values of and is ignorant 
of the spatial distribution and shapes of individual thermalized or inflating domains. One 
needs to use other considerations, specific to a particular inflationary model, to select the 
correct solution. 
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Numerical simulations of chaotic inflationary models suggested that thermalized re- 
gions always merge. However, because of the exponentially growing scales, it is difficult 
to accurately simulate the global structure of an eternally inflating spacetime, especially 
with realistic model parameters. An unambiguous resolution of the question of merging in 
a general scalar-fleld inflationary model requires further study. 
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